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ABSTRACT 

The  feasibility  of  using  x-ray  diffraction  methods  to  measure  re- 
sidual stresses  in  uranium  and  zirconium  (Zircalloy-2)  was  investigated. 
A  precision  method  was  developed  for  the  determination  of  diffraction 
peak  positions  and  the  precision  associated  therewith.   The  statistical 
tables  of  Fisher  and  Yates  were  used  to  determine  what  order  polynomial 
provided  the  best  least  squares  fit  within  the  known  precision  of  the 
observed  data.,   It  was  found  that  a  second  order  polynomial  provided  an 
adequate  regression.   With  the  aid  of  a  desk  calculator,  less  than  five 
minutes1  calculation  time  is  required  to  determine  the  peak  position  to 
a  precision  of  +  0.01°. 

A  General  Electric  XRD-3  Diffraction  Unit  was  used  in  this  investi- 
gation.  The  standard  two-exposure  method  was  employed  in  all  stress 
measurements  with  the  oblique  position  taken  at  45° .   Cobalt,  copper, 
and  chromium  radiations  were  investigated  in  determining  the  best  high- 
angle  lines  to  be  employed.   The  stress  constants  were  determined  using 
annealed  specimens  of  uranium  and  Zircalloy-2  stressed  in  bending.   The 
stress  was  measured  with  two  strain  gages  placed  one  on  either  side  of 
the  irradiated  area. 

The  stress  constant  for  uranium  was  determined  to  be  1308  +  110 
psi/0.01°  shift  in  A  2 0  for  copper  radiation  on  the  (116) planes  at 
2  9=  .158.3° o   The  stress  constant  for  Zircalloy-2  was  determined  to  be 
430  +  1  psi/0.01°  shift  inA2  9  for  chromium  radiation  on  thefl0o4) 
planes  at  29=  156.4°. 

The  surface  stress  distribution  adjacent  to  the  butt  weld  of  two 
flat  plates  of  Zircalloy-2  was  measured,,   It  was  found  that  the  weld 
produced  a  significant  effect  upon  the  stress  distribution  adjacent  to 
the  weld. 

It  is  concluded  that  the  expected  precision  as  determined  by 
statistical  analysis  can  be  attained  using  the  method  developed  herein. 
It  is  recommended  that  a  round  robin  specimen  be  sent  to  several  stress 
analysis  laboratories  to  fully  verify  this  fact  and  to  establish  confi- 
dence in  this  method  of  residual  stress  measurement. 
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L  INTRODUCTION 
Since  the  work  of  Lester  and  Aborn  (1)  in  1926,  the  use  of  x-ray 
diffraction  techniques  to  determine  residual  stresses  has  been  consider- 
ed feasible  —  at  least  from  the  theoretical  point  of  view0   Bragg 's  law 
as  pertaining  to  diffraction  of  x-rays  from  a  crystal  may  be  stated  as  j 

nA=  2  d  sin  0  (1) 

where       n  is  an  integer 

/\   is  the  wavelength  of  the  radiation 

d  is  the  spacing  between  atomic  planes  in  the  crystal 

0  is  the  angle  of  incidence  between  the  x-ray  beam  and 
the  atomic  plane „ 

For  plane  stress  as  in  the  surface  of  a  stressed  specimen,  classi- 
cal theory  of  elasticity  predicts  that? 

o-  =  4£  _E_   1 

d  1+u  sin2V  {2) 

where       £d  is  the  difference  in  the  lattice  spacing  of  the  (hkl) 
planes  parallel  to  the  surface  and  at  an  angle  V  to 
the  surface  normal 

o~  is  the  stress 

E  is  the  modulus  of  elasticity 

u  is  Poisson's  ratio 

Combining  the  above  two  equations g 

°-  =(—)    :    ffi —  -  <K>  (AG)  (3) 

V 1+u/ sin2  Y   tan  0 

where  K  is  the  stress  constant » 

1     For  complete  derivation  see  (2) „ 


Equation  (3)  indicates  that  the  stress  in  the  surface  of  a  speci- 
men is  directly  proportional  to  the  difference  between  the  peak  position 
values  of  9   as  observed  in  the  normal  and  oblique  observations .   The 
direction  of  this  stress  is  defined  by  the  intersection  of  the  plane 
containing  the  x-ray  beam  in  the  two  observations  with  the  surface  of 
the  material o 

A  great  deal  of  work  has  been  done  in  the  past  few  years  investi- 
gating mechanical  failure  of  materials  —  particularly  in  attempting  to 
correlate  the  material  strength  and  fatigue  endurance  to  the  physical 
and  metallurgical  characteristics  of  a  material <,   Most  investigators 
(3), (4), (5)  believe  the  residual  surface  stresses  to  be  an  important 
parameter  in  this  regard „ 

Several  methods  have  been  devised  to  measure  residual  stresses . 
Of  these,  only  the  x-ray  diffraction  method  is  non-destructive 0   Another 
advantage  is  that  only  the  stresses  in  the  very  surface  fibres  are 
measured  as  compared  to  the  integrating  effect  of  the  various  mechanical 
dissection  techniques  (6) .   However^  the  x-ray  diffraction  method  has 
never  become  a  generally  accepted  method  of  stress  determination  (6)  be- 
cause of  the  lack  of  precision  in  the  reported  results „ 

It  is  the  firm  belief  of  the  authors  that  the  method  itself  is 
quite  reliable o  (  It  is  rather  the  improper  experimental  approach  to  the 
problem  by  most  investigators  which  has  caused  the  difficulty „   Since 
precision  of  the  order  of  0o01°  is  desired  in  measuring  diffraction  peak 


positions,  one  must  be  extremely  careful  in  the  experimental  procedure 
and,  more  important,  in  the  treatment  of  the  observed  data.,   Precision 
is  a  mathematically  defined  concept  (7)  and  one  often  overlooked  or 
merely  estimated  roughly  in  research  as  well  as  in  engineering  work. 
With  the  exception  of  Ogilvie's  method  of  least  squares  determination 
of  the  diffraction  peak  position  (8) ,  there  is  no  evidence  of  an 
approach  sufficiently  careful  to  warrant  the  expected  and  assumed  pre- 
cision o   It  is  little  wonder  that  there  has  been  a  lack  of  precision 
and  accuracy  in  the  results  reported  using  x-ray  diffraction  methods » 
When  dealing  with  statistical  data,  one  must  apply  the  well  es- 
tablished principles  of  statistics  in  order  to  obtain  the  maximum 
correct  information  from  the  observed  data.   The  observed  counting 
rate  of  diffracted  x-rays  is  definitely  statistical  in  nature  as  was 
verified  by  application  of  the  Chi-squared  test  to  a  set  of  observed 
data.   In  experimental  design  it  is  appropriate  to  have  the  statistical 
precision  limitations  of  the  same  order  of  magnitude  or  less  than  the 
inherent  errors  in  the  apparatus 0 


lie  PROCEDURE 
A  precision  method  of  diffraction  peak  position  determination  was 
developed.   The  regression  method  by  orthagonal  polynomials  (9) ,  (10) , 
(11)  was  applied  to  observed  intensity  data  across  the  diffraction  peaks 
chosen  for  subsequent  stress  evaluation »   65,536  counts  were  observed 
at  each  of  seven  equally  spaced  values  of  angle  across  the  peako   Use  of 
the  Fisher  and  Yates  tables  (11)  provides  a  quick  method  for  determin- 
ing the  order  of  polynomial  which  should  be  fitted  to  best  approximate 

3 
the  observed  data0    The  criteria  for  fitting  is  the  least  squares  fit 

which  results  in  a  mean  square  deviation  of  the  same  value  as  the 

variance  of  the  observed  data.,   It  was  found  that  a  second  order  polynom- 

4 
ial  was  the  appropriate  function  to  be  usedo    The  direct  method  of  peak 

determination  as  recommended  by  Ogilvie  (8)  was  developed  in  terms  of 

the  orthagonal  polynomials »   This  method  in  conjunction  with  the  Fisher 

and  Yates  tables  greatly  simplifies  this  calculation  for  any  number  of 

points  (up  to  75  points)  »   Once  the  peak  has  been  ""counted,  M  less  than 

five  minutes'  calculation  time  is  required  by  this  method  to  determine 

the  diffraction  peak  position  from  seven  points o   This  least-squares 

method  was  used  throughout  this  worko 

Uranium  and  Zircalloy-2  specimens  were  irradiated  with  x=radiation 

from  cobalt,  copper,  and  chromium  in  order  to  determine  the  appropriate 

2  Section  c  of  Appendix » 

3  Section  C  of  AppendiXo 

4  Section  C  of  AppendiXo 
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5 
diffraction  lines  to  be  observed^   As  high  a  26  value   as  possible  was 

chosen  to  get  good  precision.*   The  reason  for  this  may  be  seen  as 

follows? 

From  Equation  (3) ,  for  constant  stress  and  V 

A   9  a  tan  9  (4) 

which  shows  that  a  high- angle  line  is  most  desirable  for  best  precision » 

The  lines  chosen  weres 

Uranium       Copper  Radiation   29=  158 ° 3°  (116)  planes  o 

Zircalloy-2   Chromium  Radiation  2  9=  156040  (l0o4)  planeso 

A  General  Electric  XRD-3  diffraction  unit  was  used  throughout  this 

investigation 0   The  slit  arrangement,  filtering,  high  voltage  and  current 

values  were  varied  to  determine  the  optimum  arrangement «>   In  general  it 

was  desired  to  obtain  as  high  a  counting  rate  as  possible  to  yield  ade- 

quate  precision  in  a  reasonable  timec   No  attempt  was  made  to  determine 

the  background  counting  rate  since  no  variation  was  to  be  ejected  over 

the  region  considered  (0o30°)o   The  final  experimental  arrangement 

chosen  wass 

Filter  ■==■  none 

Divergence  Slit  <==  3°  with  Solar  baffles  irradiating  an 

area  6  mm  x  1  cm  on  the  specimens „ 

Receiving  Slit  --  0ol°  (positioned  on  the  focus  circle) 

High  Voltage  —  47 „ 5  kv 


5  2 9  the  more  commonly  used  parameter  will  be  used  hereafter, 

6  See  (8)  for  complete  discussion  of  the  focusing  conditions 


Current  --  16  ma 

Oblique  Angle  —  V=  45° 

Collimating  Receiving  Slit  <==  0o3° 
It  was  found  necessary  to  allow  an  hour's  warm-up  time  for  the 
machine  with  the  x-ray  tube  on  before  a  stable  counting  rate  could  be 
observed.,^  This  is  most  important  for  precise  worko   An  appreciable 
shift  in  the  apparent  peak  position  will  result  if  the  average  inten- 
sity is  changing  during  the  period  of  time  required  to  ""count011  a  peak* 
Also,  the  diffraction  unit  must  be  carefully  aligned  prior  to  precise 
stress  evaluations o 

The  stress  calibration  constants  were  determined  by  observing 
the  change  in  ^26  as  the  surface  stress  was  varied  known  amounts „ 
Specimens  of  uranium  and  Zircalloy-2  were  vacuum  annealed  at  500°  C 
for  8  hours  to  remove  any  residual  stresses  caused  by  fabrication D 
Type  A-7  SR4  strain  gages  were  mounted  on  either  side  of  the  area  to 
be  irradiated .   Each  gage  was  read  before  and  after  each  run5  and  the 
stress  was  calculated  from  the  average  strain  value  observed  as 
measured  with  a  type  L  Baldwin  Strain  Indicator  ser  J53470o   A  Sanborn 
Twin  Visa  continuous  strain  recording  system  was  used  to  verify  that 
the  strain  remained  constant  throughout  a  run. 

The  specimens  were  stressed  in  a  four<=point  bending  device, 
thereby  assuring  a  constant  stress  over  the  irradiated  area»   This  was 
7     Verified  by  Chi-squared  testo 


verified  experimentally  by  the  strain  gage  measurements  <,   The  experi- 
mental setup  is  shown  in  Figure  1^  which  shows  the  specimen  in  place 
and  mounted  on  the  goniometer  of  the  diffraction  unite   The  alignment 
pedestal  is  shown  in  position  for  checking  alignment  of  the  specimen 
on  the  axis  of  the  goniometer „   Alignment  of  the  specimen  was  checked 
before  each  run  to  within  +  Qo001  inches0 

The  data  for  each  peak  was  taken  as  follows „   A  slow  scan  of  the 
peak  was  madeD   The  apparent  maximum  point  of  the  peak  was  determined 
to  the  nearest  0=05° „   The  time  for  16 9 3 84  counts  was  recorded  for  each 
of  seven  points  0o05°  apart  centered  about  the  estimated  maximum  value „ 
Four  separate  counts  across  the  peak  were  performed „   This  time  sequence 
is  more  desirable  than  merely  counting  each  point  four  times  and  then 
shifting  to  the  next  2  e  value,  as  any  transient  effects  tend  to  be 
eliminated,, 

The  geometrical  correction  factor  due  to  any  slight  misalignment 
was  taken  to  be  the  apparent  A2  e  value  for  a  standard  silicon  specimen 
and  a  powdered  tungsten  specimen  at  approximately  the  same  2©  angle  as 
that  used  in  the  stress  constant  determination °   Although  this  correction 
is  relatively  unimportant  in  evaluating  the  stress  constant^  K$    it  is 
of  prime  importance  in  absolute  stress  measurement «,   It  was  found 
appropriate  to  reevaluate  this  constant  for  each  experimental  setup 0 
In  order  to  verify  the  precision  method  developed^  each  author  evaluated 
the  stress  constant  on  a  different  specimen „ 
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FIGURE  I 
Calibration  Specimen  Mounted 
on  the  XRD-3  Goniometer 
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FIGURE  I 
Calibration  Specimen  Mounted 
on  the  XRD-3  Goniometer 


The  weld-induced  stress  distribution  adjacent  to  the  tungsten 
electrode  gas-shielded  arc  weld  of  two  flat  plates  of  Zircalloy-2  was 
measured*   The  1°  divergence  slit  was  used  in  this  stress  analysis,, 
limiting  the  area  examined  to  3  mm  x  1  cm  with  the  1  cm  side  parallel 
to  the  weldo   The  welded  specimen  was  mounted  as  shown  in  Figure  II 
for  stress  evaluation.   A  stress  pattern  was  determined  as  a  function 
of  distance  from  the  center  of  the  weld. 
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FIGURE  II 
Welded  Plate  Mounted  for  Stress  Analysis 
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FIGURE  II 


Welded  Plate  Mounted  for  Stress  Analysis 
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IIIo  RESULTS 


The  results  for  uranium  are  plotted  in  Figure  III,   The  stress 


calibration  constant  was  evaluated  as  1308  +  110  psi/0o01   shift  in 

Q 

A  2  6o   This  value  was  determined  by  a  weighted  least-squares  corre- 
lation with  the  predicted  linear  variation „   The  theoretical  value 

9 
was  calculated  to  be  781  psi/0o01°  shift  of  A  2  6b    The  precision  limits 

quoted  for  the  stress  constant  as  well  as  those  indicated  on  the  plotted 

points  in  the  A2  9  direction  were  calculated  from  the  known  precision  of 

the  observed  data»   The  precision  limits  for  stress  were  estimated  from 

the  +  7o5  /tin/in  observed  precision  in  reading  strain „   The  standard 

10 
law  of  propagation  of  errors  was  employed  in  these  calculations «, 

The  results  for  Zircalloy=2  are  plotted  in  Figure  IV  „   The  stress 

calibration  constant  was  evaluated  as  430  +  1  psi/0<,01°  shift  in  A20 

11 
using  the  same  methods 0  The  theoretical  value  was  calculated  to  be 

12 
377  psi/0o01°  shift  A2  9  . 

Figure  V  shows  the  variation  of  stress  in  the  welded  plate „ 

These  stresses  were  measured  in  a  direction  normal  to  the  direction  of 

the  weld  along  the  center line  of  the  specimen „ 


8  E  =  30xl06  psio 

9  u  =  0o28o 

10  Section  B,  Appendix » 

11  E  =  13o8xl0°  psi0  12     u  =  0o32 
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FIGURE  V 
Plot  of  Welded  Plate  Surface  Stress  Distribution 
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IV c  DISCUSSION  OF  RESULTS 

As  shown  in  Figures  III  and  IV,  the  results  correlated  quite  well 
with  the  predicted  linearity  of  stress  as  a  function  of  A2  e  0      In  each 
case  the  plotted  curve  does  not  pass  exactly  through  the  origin „   The 
intercept  of  the  curve  with  the  ordinate  axis  indicates  the  residual 
stress  in  the  specimen •   This  is  to  be  expected,  since  it  is  extremely 
difficult  to  produce  completely  stress-free  specimens „   The  deviation 
of  the  theoretical  and  experimentally  determined  values  for  the  stress 
constants  are  probably  due  to  the  anisotropy  of  the  material  constants a 

The  correlation  between  the  observed  and  calculated  precision  is 
excellent o   The  data  points,*  with  but  one  exception,  lie  well  within 
two  calculated  standard  deviations  of  the  least  squares  fitted  straight 
lines . 

The  precision  limits  of  stress  for  uranium  are  greater  than  those 
for  Zircalloy-2  because  of  the  difference  in  the  modulus  of  elasticity <, 
The  variation  in  the  precision  of  A2  0  is  due  principally  to  peak  broad- 
ening o   This  effect  is  included  in  the  denominator  term  of  the  error 
equation . 

It  is  interesting  to  note  that  even  with  the  greater  spread  of 
data  for  the  uranium-B  specimen,  the  precision  limits  of  the  resultant 
stress  constant  were  less  than  for  the  uraniura-A  specimen 0   Although 

13    Error  equation  is  derived  in  Section  C  of  the  AppendiXo 
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this  result  is  not  to  be  expected  in  every  case,  the  additional  infor- 
mation provided  by  the  extra  point  for  the  uranium-B  specimen  outweighed 
the  comparative  lack  of  precision  for  that  specimen o 

The  observed  surface  stress  distribution  in  the  welded  plate 
(Figure  V)  indicates  that  the  plate  was  not  fully  annealed  prior  to 
welding.   This  is  substantiated  by  two  factors e   First *  broadened 
diffraction  peaks  were  observed  in  the  regions  away  from  the  weld5  the 
condition  normally  present  in  an  unannealed  plastically  deformed  material  <, 
Second,  large  compression  stresses ,  which  are  characteristic  of  rolled 
plate,  were  observed  far  removed  from  the  welde   The  observed  stress 
pattern  is  presumed  to  be  due  to  the  combined  effects  of  fabrication  and 
welding o 
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Vo  CONCLUSIONS  AND  RECOMMENDATIONS 
It  is  concluded  that  with  the  proper  experimental  approach  and 
statistical  analysis  of  the  observed  data,  the  x-ray  diffraction  method 
for  measuring  residual  stresses  is  indeed  a  practicable  oneQ   This  fact 
is  especially  emphasized  by  the  excellent  results  herein  obtained  using 
diffraction  peaks  that  would  hardly  be  considered  sharp  peaks,  as  shown 
in  Figure  VI.   The  excellent  correlation  between  the  results  from  the 
A  and  B  specimens  confirms  the  validity  of  the  method . 

It  is  recommended  that  this  fact  be  firmly  established  by  start- 
ing a  round-robin  specimen  to  various  research  and  testing  laboratories . 
This  specimen  should  be  accompanied  by  instructions  for  the  method  re- 
commended herein  in  determining  the  residual  stresses  in  the  specimen . 
An  alternate  procedure  would  be  to  re-evaluate  the  data  and/or  re-measure 
the  stresses  in  previously  circulated  specimens  (6) 0      In  this  way  the 
true  merit  of  the  x-ray  diffraction  method  would  soon  become  apparent o 
It  is  also  recommended  that  further  investigation  be  carried  out 
in  the  fields  ofs  — 

lo  a  comparative  investigation  of  higher  intensity  lower 
angle  lines o 

2o  a  comparative  investigation  of  photographic  techniques 0 

3o  a  complete  stress  analysis  of  the  weld-induced  stress 
in  a  properly  annealed  specimen., 

There  are  numerous  possible  applications  of  this  technique  in 
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the  field  of  reactor  design „   Some  of  these  would  includes 

lo  Evaluation  of  stresses  produced  by  local  hot  spots  in 
Zircalloy=»2  cladding „ 

2.  Determination  of  the  stresses  induced  by  welding  and 
fabrication  techniques  in  fuel  elements „ 

3  o  Determination  of  the  appropriate  heat  treatments  to 
properly  relieve  residual  stresses ° 

The  curve-fitting  and  error  propagation  techniques  described 
herein  are  perfectly  general 0   They  can  and  should  be  applied  in  any 
experimental  investigation,  in  order  to  properly  establish  the 
precision  of  the  results . 
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VI o    APPENDIX 
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REGRESSION  AND  ERROR  ANALYSIS 
A.  Regression  Analysis 

Regression  is  the  determination  of  the  coefficients  of  an  assumed 
mathematical  expression  purporting  to  represent  the  true  means  of  an 
array  of  data.   The  selection  of  the  regression  equation  may  be  left  to 
the  judgment  of  the  analyst,  in  which  case  the  simplest  form  producing 
adequate  results  is  desired,  or  there  may  be  physical  laws  which  will 
govern  the  selection.,   Whatever  form  the  regression  takes  it  must  be 
remembered  that  the  methods  presented  here  are  only  mathematical  means 
to  an  end. 

Historically,  the  name,  regression,  originated  with  Galton,  who 
first  applied  the  method  to  determine  how  a  population  could  remain  in 
dynamic  equilibrium  if  the  sons  inherited  the  characteristics  of  the 
father.   He  found  that  the  sons'  deviation  from  the  mean  stature  was 
less  than  the  fathers'  and  consequently  the  sons  ""regressed00  toward  the 
mean. 

Consider  an  array  of  data  with  n  columns  and  m  rows  such  that  one 
of  the  variables  may  be  taken  as  independent.   The  first  and  probably 
most  useful  relationship  between  the  variables  would  be  that  they  were 
linearly  related  or  nearly  so.   The  questions  to  be  answered  ares  that  if 
they  are  not  exactly  linearly  related,  what  degree  of  approximation  must 
be  made  to  obtain  such  a  relationship?,  and,  is  the  difference  from 
linearity  significant?   Assuming  X  to  be  the  independent  variable,  and 
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Y  the  dependent  variable,  the  coefficients  of  the  equation 

y'  =  ax  +  b  (5) 

are  to  be  determined .   y  .  is  the  variate  in  the  i  th  column  and  the 

j  th  row.   The  average  value  of  the  i  th  row  is  y  , 

i 

j=m 

Y.  -1  I       y.   .  (6) 

For  reasons  which  will  become  obvious  it  is  easier  to  determine  the 

coefficients  of 

y'  =  a  (x  -  x)  +  b  ,  (?) 

where  x  in  the  usual  notation 

n 

x  =  1      i 

n     '  (8) 

An  appeal  is  now  made  to  the  Principle  of  Least  Squares  upon  the  assump- 
tion that  it  will  lead  to  the  most  probable  values  of  the  coefficients „ 
Therefore  the  minimum  of 

n  2 

G  (a,b)  =  X  lY±   ~   y'J  (9) 

1    x  1 

is  desiredo   Substitution  of  the  assumed  function  yields 


n 
(a,b)  =  Y.    [y±  -  b  -  a  (x±  -  x)J 


:io: 


The  minimizing  conditions  are 

n 


iG  = 


-  -2  Z  [y±  -  b  -  a  (x±  -  x)l  =  0 


11] 
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6G 

da 

which   leads   to 


n 
-2  I 


1  L 


a  (x.  -  x) 
i 


(xi  -  x) 


n  n 

2  y .  =  bn  +  a  J  (x 
1  1 


x) 


n 


n 


n 


5"  (x,  -  x)  y.  -  b  T  (x.  -  x)  +  a  7  (x 


-  o  , 


x) 


(12: 


13' 


(14) 


Since  by  definition 


n 


and 


I 

(x.    -   x) 

=   0 

1 

l 
n 

b 

_  i 

n 
n 

=  y 

I    (x 
1 

-   x)    y. 

l 

a 

5  2 

I  (x.  -  x) 

1 
with  the  resulting  desired  expression 


y  =  (x 


n 

I 

*\    1  (x- 


x)  y 


£  i, 


n  2 

I    (x.  -  x) 


15! 


:i6) 


which  is  defined  as  the  regression  of  y  on  xe 

With  the  constants  of  the  equation  now  determined  the  question  of 
the  justification  of  a  linear  regression  may  be  undertaken =   The  usual 
approach  is  to  separate  the  deviation  of  any  variate  from  the  mean  into 


-  26  - 
two  parts,  the  deviation  due  to  regression  and  the  deviation  from  the 
regression.   There  follows 

y  -  y±     =  y  _  y^  +  y^  -  y± 

n         2    n  r  1  2 


I  (y  -  y.)  =  I   [(y  -  v[)  +  (y[  -  vj\ 


n  2     n                     n  2 

=  I  (y!  -  y)   +  2  J  (y|  -  y)  (y.  -  y!)  +1  (y  -  y! ) 

1    1  1    1        i    !    x  i    i 

n         2  n         2   n          2 

I  (y  -  y.)  =1  (yf  -  y)  +  I  (y.  -  y1)  Mft% 

1       i  t   i         i   i    i  (18) 


since 


n 

X  (y!  -  y)  (y  -  y8)  =  o.  (19) 

1   i        i    i 

Thus  the  important  result  that  the  variance  of  the  mean  can  be  summed 
from  the  variance  due  to  regression  and  the  variance  from  the  regression0* 
The  variance  from  regression,  (residual) ,  is  certainly  a  good  measure  of 
the  justification  of  the  linear  regression  as  compared  with  some  other 
assumed  mathematical  expression „   Although  there  are  no  set  rules  it  would 
seem  inappropriate  to  select  an  expression  of  higher  order  if  the  average 
variance  from  regression  is  appreciably  less  than  the  average  variance 


*  It  is  entirely  likely  at  this  stage  that  an  estimate  of  the  statistical 
significance  of  the  regression  is  desiredo   This  may  be  obtained  since 
there  are  now  available  two  separate  estimates  of  the  same  variance  with 
n^L  and  n   degrees  of  freedom.   The  "Z"    or  its  equivalent  "F"  test  or  the 
better  known  Student's  "t10  test  may  be  appliedo 
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of  the  means  of  the  variates.   Mathematically  this  may  be  expressed  ass 

1  2      1     2 

n  2.   i    i   -  «  2-      yi  (20) 

n  n 

2 
where  6     is  the  variance  of  the  mean  of  the  i  th  column  in  the  array „ 
yi 

If  linear  regression  is  found  to  be  unsuitable 3    a  logical  exten- 
sion is  to  investigate  polynomial  regression.   A  straightforward  approach 

would  be  to  assume  an  equation  of  the  type 

2  n 

y  =  b  +bx  +  bx  +...bx 
o    1     2  n 

and  proceed  as  before  by  taking  partial  derivatives  and  applying  the 
minimizing  conditions,   A  better  approach,  as  will  be  seen,  is  to  define 

a  set  of  orthagonal  polynomials 

2      3      n 

Z  =  a  +  a,x  +  a  x  +  a  x  . +.  a  x  12\\ 

n        1     2      3      n 

where 


n 

2 

n 


I  z„  =  0  i  (22; 

1 


n 

I  \    =  °  *  (23] 

1  K 

n 

£  Z  Z  =  0,  n  ft   k0  (24] 

1  k  n 


By  this  definition 


(25; 
(26; 


28 


Z0  =  (x3  -  3n  -7  x) 
3  20 

2 


=  [x4  -  3n  -13  x2  +  3(n  -1)  (n2-9) 
L       14  560 


(27; 


28] 


4    L       14  560 

As  an  example  consider  fitting  a  second  degree  polynomial  to  the 
data.   The  assumed  expression  is 


y'=A  +A  Z  +  A  Z. 
°     1   1     2   2 

The  function  to  be  minimized  isi 


;29] 


£  2 

G  (Ao'  V  V  "  2  (yi  "  y±} 


Upon  substitution 


G  (A  ,  A 
o 


and  differentiation 


1 
n 


.*    A   )    -  J    (y      -   A      -   A  Z 

12  l         i  °  1    1 


-    A   Z    )       , 
2    2 


JG_      =   2 
dA_ 


n  r 


A      -    A   Z      -   A   Z    )  (-1) 
°  11  2    2 


]■• 


(30) 


:3ii 


(32: 


n 

dA.  f 


(y.    -   A      -A   Z      -   A   Z    ) (-Z    )     1=   0 
Li  oil  22  1J 


6g 


n    „ 
=   2£ 


-2  1    L      x  o  1 

By  equations     (22),     (23),    and    (24) 


(y     =-A     -   A  Z      ~AZ)  (-Z    ) 
i  o  i    1  2    2  2 


)]  =   0 


34: 


n 

X  (y    -  a) 

1         i 
n 

n 


=   0 


=  y 


;3s: 


:36! 
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n 

2 

zi  ■ 

n 

\ 

I 

I 

1 

n 

1 

Al 

Iz 

1      1 

y± 

n 

? 

37) 
Z  Z^ 
1   1 


n       n 

^12     1   Z   i 

n 

I  Z  y 

A  =  1~  2  >  (38) 


2     n   2 

2    I  z2 

1 

Further  generalization  leads  to 


n 

1   n  i 

i 

n 

n       2 
XZn 

1 

(39; 


As  mentioned  previously  the  advantage  of  the  method  can  be  readily 
seen,  since  the  sums  of  the  squares  from  regression  ares 

Linear    §  AlZy  (40) 

Quadratics    AXZ  y     +AIZy  C41) 

1      1   i  2      2    i 

nth  sAZZy      +  A  2  Z  y      +    *  0  o    +  A  XZ  y  (42) 

lli  2      2i  nnJi 

The  sum  of  the  squares  from  regression  can  be  easily  determined  and  the 
effect  of  each  additional  stage  of  fitting  evaluated,. 

There  still  remains ,  however,  the  determination  of  the  constants 
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of  the  orthagonal  polynomials.   The  advantage  of  using  orthagonal  poly- 
nomials in  the  general  case,,  when  unequally  spaced  and/or  unequally 
weighted  values  of  the  variates  exist,  is  difficult  to  justify  in  so  far 
as  numerical  labor  is  concerned^  but  for  equally  spaced  and  weighted 
values,  the  coefficients  of  the  polynomials  may  be  determined  once  and 
for  all,  and  tabulated,  for  a  cardinal  number  spacing .   Tables  of  the 
polynomials  for  orders  up  to  and  including  the  fifth  degree  for  up  to 
and  including  seventy-five  points  have  been  developed  by  Fisher  and 
Yates o   These  tables  also  include  values  for  computing  the  sums  of 
squares  necessary,, 

B.  Precision  and  Propagation  of  Errors 

Consider  the  problem  of  calculating  the  precision,  or  error,  of  a 
function  made  up  of  independently  measured  functions  (i0e0,  empirical 
data)  whose  errors  are  known »      Such  a  problem  would  be  the  calculation 
of  the  total  resistance  coefficient  for  a  model,  where  the  resistance, 
speed,  and  surface  area  cannot  be  exactly  determined 0 

The  required  assumptions  ares 

a=   Each  quantity  to  be  considered  must  be  arrived 
at  independent ly . 

bo   The  errors  in  the  measured  quantities  must  be 
small o 

Co   Positive  and  negative  errors  must  be  equally 
probable o 

Let  the  true  values  of  the  measured  quantities  be  x. ,  x  8    x  o  etc. 

1   2    3 
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These  values  are  to  be  combined  in  some  manner  to  give  the  desired 
function 

X  =  f  (xx,  x2,  x3,  .  ..)  .  (43) 

Let  the  true  errors,  which  are  unknown,  of  the  measured  quantities,  be 
e^,  e_,  e3  „..  etcc   The  unknown  error  in  X  will  be  E.   Then 

X  +  E  =  f  (x,  +  e  ,  x  +e,x  +e   <,<>.)  <,  (44) 

By  the  required  assumption  (b)  e,/x.    is  much  less  than  one,  and  a 
Taylor's  Series  Expansion  of  equation  (44)  may  be  made  with  only  the 
first  order  terms  retained,, 

X  +  E-X  =  E=  (dx/dx  )  e  +  (dX/dx  )  e  aooo  (45) 

By  squaring  both  sides  of  equation  (45) 

2  2   2  2   2 

E   =  (dx/dx  )   e,   +  (dX/dx  )   e   +  „ . . 
11  2     2 

+  (dx/dx  )  (ox/ax  )  e  e   000  (46) 

1        2    12         2 

where  the  terms  contain  the  errors  in  the  form  e.   and  e.e  0      Since  (by 

1       i  k 

assumption  c)  the  errors  are  equally  probable  in  a  positive  and  negative 
sense,  the  cross  product  terms  will  tend  to  cancel  out»   Therefore,  the 
desired  result,  which  is  commonly  termed  the  Law  of  Propagation  of  Errors, 

is, 

2  2   2  2    2 

E   =  (dx/dX;L)   e    +  (dx/dx  )   e    +  ,OOOOB  (47| 

The  error  e  of  a  value  x  is  usually  encountered  in  the  form 

x  +  e 

This  implies  that  there  is  a  50%  probability  that  the  true  value  of  x 

lies  between  x  plus  e  and  x  minus  e„ 
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C,   Example 

Consider,  as  an  example^  the  intensity  data  from  an  x-ray  diffra= 
ct ion  peak  from  a  Zircalloy-2  specimen  under  stresso 


TABLE  I 


X-ray  diffraction  data 


■i              \ 

1 

2 

3 

4 

5 

6 

7 

J-             / 

/ 

*>  C\           ^ 

156. 50** 

156.55° 

156.60° 

156.65° 

156.70° 

156.75° 

156.80° 

J 

Z   U           / 

Time  for 
16,384 
counts 

(in  sec) 

36.22 

35.68 

34.57 

34.75 

34.94 

35.58 

36.40 

1 

36.51 

35.50 

34.25 

34.74 

35.23 

35.20 

36.69 

2 

36.56 

35.42 

34.87 

34.27 

34.83 

36.22 

36.49 

3 

36.91 

35.64 

35.34 

34.69 

34.36 

35.71 

36.26 

4 

V          \ 

36.55 

35.56 

35.01 

34.61 

34.84 

35.68 

yi    ; 

36  .46 

In  this  array 


m  =  4 

n  =  7 


Taking  Z»  values  (Z«  =  h  Z)  from  the  Fisher  and  Yates  tables  for  n  -  7 
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TABLE  II 


Coefficients  of  the  orthagonal  polynomials  for  n  =  7 


i 

Zi 

Z2 

Z8 
3 

Z« 

4 

Z» 
5 

1 

-=3 

5 

=1 

3 

-1 

2 

-2 

0 

1 

-7 

4 

3 

-1 

=  3 

1 

1 

-5 

4 

0 

-4 

0 

6 

0 

5 

1 

-3 

-1 

1 

5 

6 

2 

0 

-1 

-7 

-4 

7 

3 

5 

1 

3 

1 

I  <z')2 

28 

84 

6 

154 

84 

h 

1 

1 

1/6 

7/12 

7/20 
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TABLE    III 


Calculation  of   variance 


i 

yi 

2 
Yi 

z\*i 

Z8  y 
2  *i 

Zj  y 
3  J± 

Z«  y 
4  i 

Z»  y 
5  i 

1 

36o55 

1335.9025 

=109.65 

+182.75 

-36.55 

+109.65 

-  36.55 

2 

35.56 

1264.5136 

-  71*12 

+  0 

+35.56 

-248.92 

+142.24 

3 

35o01 

1225o7001 

-  35.01 

-105.03 

+35.01 

+  35.01 

-175.05 

4 

34.61 

1197.8521 

0 

-138.44 

0 

+207.66 

0 

5 

34o84 

1213o8256 

34.84 

-104.52 

-34.84 

+  34.84 

+174.20 

6 

35o68 

1273.0624 

71.36 

+   0 

-35.68 

-249.76 

-142.72 

7 

36.46 

1329.3316 

109.38 

+182.30 

+36.46 

+109.38 

+  36.46 

I 

248o71 

8840c 1879 

-   0.20 

+  17.06 

-0.04 

-   2.14 

-   1.42 

-    Iyi 

y  =  — —  =   35.53 
n 


y     =    12623809 
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TABLE  IV < 


Analysis  of  variance 


Deqree   of  Fit 

Sum  of   Sauares   =    (SS^ 

Residual 

Average 
Residual 

I(Y±)2    -    I(SS)i 

0 

ny2   -   8836.6663 

3 ■ 52 16 

1 

[*<zii*i>]2 

*—-  =    .0014 

I(Zil)2 

3o5202 

0o503 

2 

[I<Z2i   Vi)]2 

L             .     ,  x  J       -    3    4648 

0o0554 

O0OO8 

I(z^.)2 

3 

r^(z3i  y.>i2 

L           -3i      l    J       _      0003 

0.0551 

0.008 

I(z;.)2 

4 

[I(Z4iV] 

2 

=        n?Q7 

0.0254 

0.0036 

I.(z»  ')2 
4i 

5 

[2<z5i^]2 

0o0014 

0o0002 

I(z«   )2 

5i 
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As  discussed  above,  the  residual,  (variance  from  regression) ,  is 
the  required  measure  of  goodness  of  fitc   If  a  seventh  degree  equation 
had  been  chosen  the  residual  would  be  zero,  as  the  curve  would  pass 
exactly  through  each  point » 

However  the  precision  (ioea,  the  variance)  to  which  each  y. 
value  is  known  is  easily  calculable.,   For  counting  a  constant  number 
of  counts  (s)  in  time  (t) $    the  variance  (52)  is  given  bys 

62  =  1  (t2/s)  (48) 

m  x   ' 

(The  variance  of  the  mean  of  four  counts  is  one-fourth  the  variance  of 

any  onec)**  Therefore s 

for   m  =  4 

t  <*   35 

s  =  16384 

52  =  o019  sec2 

Thus  it  is  seen  that  the  second  order  equation  is  the  appropriate 

regression  to  be  employed. 

The  information  that  is  desired  is  not  specifically  the  regression 

equation  but  its  maximum  value  which  may  be  determined  as  follows  § 

y  =A+B9Z»+C"Z»   (in  the  notation  of  Fisher  and  Yates .] 
e  12 


**  EVANS  —  pp„  756  and  798, 


where       B 
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i 


c<  =  -1  Yi  aa 


I  (z 

2i 


^  =  0  =  B»  £?i  +  C  d22 
dx  dx        dx 


Z»  =  h  x 
1    1 


z;-h,(3c2-flif-L) 


dZ' 

dx     1 


dZ» 

— =■  =  2  x  h 

dx         2 


x 

.  .    max 


=  ^1   -^yiZli  2(Z21)2  ^_ 

C'2hn   2y.Z»    y(Z"  )2  2h 

2   *Ji2i   A  li  2 

Determination  of  the  precisions 


Assuming  y  =  y  +  e 
i    i  —  i 


x     +  E  =  x     + 
max  "       max 


[fc  '$ 


2 

_  -Ifii   h, 


:  ■  *"  ^2   vT"    ^is  independent  of  y) 
li 


Assuming  e  =  e  =6 


50! 
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E   =  6    2 


/dx  V 


dx 

dyi 


=    IT 


f^i  Hi)  zli +  (*  yizli)  z2i] 


(ly.z'J 


ay± 


i   2i 


2       .2 


3>        2 


(Iy.Z«     )       Z»       +    (Iy.Z»    )    Z' 
i   2i  ii  i    li        2i  _ 


[    Z1       Z«       =   0 
li      2i 


E   = 


ir& 


(Zy^j  ! 


ik 


(Iy,z»  )      I(z»   )     +  (ly  z»   )     i(z»   ) 

1   2i  li  i    li  2i 


E  ■ 


TTy^P 


irzI(z«    )2  +  x  Z(Z«    )2 

li  max  2i 


-  i/ 


:si] 


for  n  =   7   from  Fisher   and  Yates 


h     =  h     =   1 
1  2 


I  (Z«    )2   =   28 
li 


2(Z«    )       =84 


TT    = 


84   x   1  m   3 

28   x   1   x  2  2 


E   = 


I  y.z° 

*i    2i 


T 


9  2  \^ 

7   28   +   x   ,      84 
4  max        / 
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for  the  particular  example 


c    0-C--2)   84  1=  O0176 
max    17c06   28   2 


6  =  .137 


E  = 


137 


^  x  28  +.0003  x  84 


(17.06; 

-    .00406 

In  terms  of  spacing  which  is  .05° 

x    =  .00088  +  .000203° 
max         — 


% 


This  calculated  error  due  to  counting  statistics  is  negligible  as  compared 
to  the  geometrical  precision  of  the  diffraction  unite    The  instruction  manual 
for  the  diffraction  unit  quotes  a  precision  value  of  +  0„002°„    Therefore,  the 
final  value  of  peak  position  iss 


2  8=   156.650  +  x 


2  & 


max 


156.651  +  0,002° 


thesH817 
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